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Abstract 

We present a coupled decreasing sequence of random walks on Z (X. , i £ N) that dominates the 
edge process of oriented-bond percolation in two dimensions. Using the concept of "random walk in 
a strip " ,we construct an algorithm that generates an increasing sequence of lower bounds that con- 
verges to the critical probability of oriented-bond percolation. Numerical calculations of the first ten 
lower bounds thereby generated lead to an improved,ie higher, rigorous lower bound to this critical 
probability, viz. pc > 0.63328. Finally a computer simulation technique is presented; the use thereof 
establishes 0.64450 as a non-rigorous five-digit-precision (lower) estimate for pc- 
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1 Introduction 

Oriented percolation in two dimensions or discrete time contact process on Z is a one-parameter family 
of discrete time stochastic processes defined on {0, 1}^; the parameter p, the infection rate, taking values 
on [0, 1]. It is a well established fact that this family exhibits a fase transition, as the value of p increases 
from to 1: if p < Pc ^, the process dies out almost surely; whereas if p > pc, the process survives with 
positive probability (see P for instance ). 

As usual, in critical phenomena theory, an analytical expression for pc is unknown and its value has 
been estimated both in mathematical and physical literatures ^ Sec. 6]. 

Up to the present moment, the best (rigorous) lower and upper bounds for pc are respectively 0.6298 
and 2/3 2. 

In this paper, we present an algorithm that generates an increasing sequence of lower bounds that 
converges to the critical probability of oriented bond percolation in two dimensions and, calculating the 
first ten lower bounds thereby generated, we were able to improve the best lower bound known up to 
now from 0.6298 to 0.63328. 

'Supported by CNPq (grant N. 301637/91-1) / e-mail: belitskyQime. usp.br 
^Supported by a grant from CNPq / e-mail: thomasl@ime.usp.br 
^ "pc" is known as the critical probability of the family. 



More specifically: 

1. a numerical sequence of lower bounds for pc, {Pc'^jieN, was constructed in the following steps: 

(a) Associated to a decreasing^ sequence(in i) of random walks on Z, {X^^}net-i, that dominate 
the edge process of oriented percolation, we constructed a sequence (in i) of finite, discrete 

time Markov chains, {Yn''}n&i, whose transition probabilities are specific'^ functions of p, the 
infection probability; and 

(b) calculated the mathematical expectation of the random variable 
E (yX^^^ — Xq\Yq''^, the mean jump on the y-configuration, with respect to tt^*', the sta- 
tionary measure for F^*^ . 

(c) pf* was then defined to be the (only) value of p that nullifies the above expectation. In other 

words, when p = pc \ the random walk {X^ }neN has zero mean drift. 

2. the numerical sequence {p^c^}i^n was shown to converge in a non- decreasing fashion to Pc, 
ie Pc^ / Pc- Moreover, 

3. the first ten lower bounds {pi^\pi^\ ...,pf ^) were numerically calculated, thereby improving the 
best rigorous lower bound known up to the moment from 0.6298 to 0.63328. 

In the last section we present a simulation technique which in our opinion has important advantages in 
comparison with the usual so called Monte Carlo simulation techniques, in the sense that it exhibites a 
clear cut off between the subcritical and supercritical phases; thereby enabling a precise estimation of 
the critical probability of Oriented Bond Percolation without the aid of scaling techniques. By means 
thereof a lower bound for pc was obtained within a precision of 5 digits, viz p^^^^ = 0.64451, so that we 
can (non-rigorously) state that Pc = (0.64451 ± 0.00001). 



2 Definitions and Constructions 
2.1 The Enviroment 

Let G = {V,£) be the oriented graph, having V = {(n, m) : n G N and (n + m) is even} as its set of 
vertices/sites, and £ — {ejj„, e5^„j : (n, m) G V} as its set of bonds. 

Bond ,„ points from site (n, m) to site (n + 1, m — 1), whereas bond e!^ „^ points from site (n, m) to 
site (n + l,m + 1). Sometimes the natural association: Z <-> — 1 / r +1 will be assumed through out 
this text. 

It is useful to think of n as a (discrete) time coordinate and of m as a (discrete) space coordinate of 
the graph Q. 

Vn denotes the n-th slice of V, ie V„ = {(n, m) e V : n is fixed }, and Z„ the set of integers m such 
that (m + n) is even. 

It will be often useful in the forthcoming definitions to identify V„ as Z„, and to think of {0, 1}^" as be- 
ing {0,1}^". Therefore arises the loose, yet natural, notation: 

r]„{m) = T]n{n, m), ?7„ e {0, 1}''^". 

^in a sense to be specified later by means of coupling 
^polinomial 
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Figure 1: Graph Q, whereon the processes will be defined (Sec. 12 .31 below). 



2.2 The Probability Structure 

Let {Cnm ■ {n,m) G V,j G {Z,r}} be a family of independent and uniformly distributed 
( onto [0, 1] ) random variables, defined on the same abstract probability space f2. Starting from this 
three-indexed countable family, we construct the four-indexed uncountable family of iid Bernoulli random 
variables on the same probability space : 

{p^,™ : "i) e V, J e {/, r},p e [0, i]} (1) 

such that 

It follows straightfowardly from ( ^ that 



n p£.im = 1) = P (a) 
np^im=Q) = I-P -.^ q {b) 



(3) 



and the Fundamental Coupling: 

Pi<P2^ p^^imi^) < Vw e n, (n,m) e V, J e {l,r} (4) 

We observe that p^n^ = 1 is interpreted as a channel open to infection propagation from site {n, m) to 
site (n + l,m ~- 1); p^,'^„j = 0, as a channel obstructed to infection propagation from site (n, m) to site 
{n + l,m + 1); and so forth ... 

2.3 The Processes 

2.3.1 The Strengthened Discrete Time Contact Processes (SDTCP) 

77,1 G {0, 1}^", is interpreted as a infection state on slice V„: 

Vnifn) = 1: site m is infected at time n , , 

Vn{fn) — 0; site m is healthy at time n 
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Particulary rjQ £ {0,1}^° will denote an initial state of infection over Vo; the set of even integers, 
according to Section 12. II above. 

At this point, the following notation (to be used throughout this paper) should be kept into account: 
N = {0, 1, 2, • • • } and N = N U {oo}. 

Now for each i e N, we define (by induction on n) the sequence of {0, l}^"-valued random variables 

{ pXn }neN as 

Definition 2.1 . 

(a) pXq'' ■= Vo, so that pX^^m) — rjQ{m) ,\fm G Zq 

(b) pX^+i(to) = sup I pXi"\m - 1). pCm-i; pXL'\m + 1). p^i,,„+i; 1^ ,-m>2«}}' 
where pX^+i = sup„gz„;iG{-i,+i} {(^ + j) ■ pXi*'(m).^^,„ = l| 

The role of the indicator function in Definition |OJb) above is to infect by force all the sites lying farther 

(i) 

than 2i on the left side of pX^^^i, the utmost right infected site at time n + 1. It is natural, therefore, to 

call the stochatic process { pXn ^}n(^N, defined above, as the Strengthened Discrete Time Contact 
Process (of i-th order and infection parameter p ). 
It follows directly from Definition 12 . II that 

(a) when i — oo, the indicator function does not act anymore, and we recover the ordinary Discrete 

Time Contact Process, also called Oriented Percolation in Two Dimensions, described in 
Sec. 2]. So forth in this paper we shall refer to it as { pXn°°'^}neN ; 

(b) the sequence of SDTCPs {pX^n^jn gN is decreasing in i, in the following sense: 

i<j;i,jeN^ pXi')(m)M> pXi,'\m)[Lo],yLuen,{n,ni)eV,pe [0,1] (6) 
In particular, the Oriented Percolation Process is the weakest of them all. 
Inequality above is called Coupling of First Kind 

(c) the family of SDTCPs { pXi'^}„e N is increasing in p in the following sense 

Pi<P2^ p,Xi'\m)[Lu] < p,X^p{m)[uj],\fLj G n,{n,m) G V, i G W (7) 

Inequality |3 above is called Coupling of Second Kind 

In this paper, unless otherwise stated, 770 will be chosen to be l,<o G {0, 1}^°, where 

l.<o(m) = ( J'""^" (8) 

- ^ ' [ 0, TO > 

We shall write { "^X^'^} neN to emphasize that 770 — r], when ry is some specific element of {0, 1}^" ( 

probably different from l,<o)- Therefore the symbols {pXi'-'}„gN, { ''''-pXn^}ni£jq, ^'-pXii'\ p^i*^ have 
all the same meaning and, for simplicity's sake, the last will be chosen, when no confusion may arise. 

As usual in Particle System's notation we write rjn > 0n, whenever r]n{m) > 0„(m) for every m G Z„. 

Definition 12. II also implies that, for r],6 e {0, 1}^" 

V>9^ 1X^^ (m) [uj] > Ixi'^ (m) [lj] , Vp G [0, 1] , (n, to) G V, w G 1^ , i G N (9) 
Inequality |^is called Coupling of Third Kind. 
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Figure 2: A realization lu of tlie pXP' process . Full-black arrows are open, whereas dotted ones are closed to 
infection propagation. Black sites/cells are infected, white ones are healthy, while the gridded ones were reinfected. 



2.3.2 The Right Edge Processes (REP) 

Given a particular SDTCP pXi^'' and assuming that p^o*'' 0, the second line of Definition |OJb) 
above defines a (non-markovian) random process on Z denoted by { pX^ }neN ( or pX ^ in abbreviated 



fashion). In Figure [51 above pX^ = 0, pX^ = —3, pX 



-Cr(2) 



-2, pX 



-Ef(2) 



3, riX A 



Again, in case of « = oo, the Edge Process of Oriented Percolation cited in jlj is recovered. Throughout 
this text, this term will be employed in a generilized way (i G N). 

It is useful to think of pX^ as random walks on Z. 
2.3.3 The Induced Markov Chains 

Definition 2.2 The Markov chain { pYn^^^n&ij * G with state space 
is called Induced Markov Chain(IMC) for the SDTCP pXi''' . 

In what follows, a generic element a e iS*^*-*, will be labeled by n e {0, 1, . . . , 2' — 1} according to 



Accordingly, in Figure [3 above pY^^^ = a^, pY^''' = a-i, pY^^' = dg, pY^^' = a^, pY^^' = cti. 



(2) _ 



.(2) _ 



.(2) _ 



42) _ 



The transition probabilities for the IMC pY, 



1l 



Yif> =ai ,0<l,m<2' -1 



(10) 



(11) 
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are all strictly positive polinomial functions of p (provided < p < 1). Therefore Tr^*^ its stationary 
measure on 5*^*^ , is well defined. 
The notation 

is self explanatory. 

At this point, in accoradance with Sec.3.1], we state 



Definition 2.3 



PtLk)ip) F = ,x':^ + (1 - 2fc), = am I pFi') - m 

the transition probabilities from state ai to state am with a jump of magnitude (1 — 2k) 
Definition 2.4 

oo 2'-l 

Ml^\p) E E (1 " ^k).Pl;l^,^{p), 0<1<T-1 

the mean jump/drift of the REP on configuration ai 
and 

Definition 2.5 



A: e N,0 < < 2* - 1 



the mean jump /drift of the REP. 



M^''\p) is also called the (right)edge speed in accordance with Section l2.4l bellow. As it will soon become 
clear, Definition 12. 51 is of fundamental importance in this paper. 

Under this framework, the SDTCPs pXi^\ described above, can be regarded, in case of i G N, as 
Markov Chains/Random Walks in a 2^-rowed Strip: 



a slightly different idea of Markov Chains in a Half Strip ^ developed in [31 Sec.3.1]. 
2.4 Critical Probabilities 

In Section |21below it will be shown that, for i G N and p G [0, 1]: 

(i) M'-*-'(p) is a scrictly incresing function of p; 

(ii) M'^^\p) has only one real root ( in (0, 1]); 

(iii) ^nM«b) as 
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In case of i = oo, a{p), the (right) edge speed of oriented percolation, plays precisely the same role of 
M^^\p) in the finite case depicted above^.(For details, see P, for instance). So the notation 

suggests itself and we state the following 

Definition 2.6 The critical (infection) probability p^f' for the family of stochastic processes 
is the only real root of M^^^^{p) (in (0, 1] ), the edge speed. Hence 

M(*) = 0, Vi e N 

The heuristic meaning of Definition 12.61 above is : 

(i) for p < p^c \ lim„_>oo p^i^ = ~oo a s, ie the infection dies out with probability one; 

(ii) for p > pc \ lim„_,oo pX'^n — the infection spreads out over all Z. 

In the sequel, we prove an important relation concerning the pi*^s, i € N, just defined, viz. pc''' i Pc°°^ ^• 
This non-decreasing convergence to the critical probability of Oriented Percolation (to be called The 
Convergence Theorem), besides the possibility of 
Section 0] below), is the cornerstone of this work. 



Convergence Theorem), besides the possibility of calculating the pi^\ (i G N) by algebraical means 



3 The Convergence Theorem and PreUminary Results 

Lemma 3.1 ^-^-^ M^^ip) a s, Vi G N, p G (0, 1] 
Proof: 

First case (i £ N): Let nf^ be the (random) number of visits that the IMC pY^^^ makes to state u; up to 
time n ( so that X]f=o^ '^i*' ~ '^) ^^'^ Ji^l the fc"* jump of the REP pXi'"' on state/row oi. 
The (strong) Markov property of SDTCP pXl'' makes : k £ {1, 2, 3, . . .} iid rvs with E [j/^'j = m/''. 
Now 

p^n _ \ ^ \ ^ ■^i.fc _ \ ^ \ ^ "^i.fc _ \ ^ \ ^ ■^i.fe 

n ~ 2^ 2^ n ~ ^ ^ ri^^) ' n ^ ^ n ^ „(') ' 
1=0 fc=i 1=0 fe=i 1=0 fc=i 'n 

!h_ n vr« as and y 4^ ^ M^^Hp) a s 

Taking the limit n ^ oo and bearing definition 12.51 in mind, we get the desired result. 
Second case (i — oo): See ^ pag.1004] 

□ 



= — oo,when p < pc, so that the strict increasing behaviour does not apply to ol(p) precisely. 
^Pc""' Pc, putting the notations of Sections 1 and 2 into agreement. 
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Lemma 3.2 For i £ N , each function M^*-' : (0, 1] — > (— oo, 1] is strictly increasing (in p). Moreover, 

p M(')(p) 

is a surjection from (0, 1] onto (— oo, 1]. 
Idea of Proof: 

The non-decreasing behaviour of M^''\p) follows from the Second Kind of Coupling (inequality |7|above) and 
Lemma 13.11 iust established. Definitions 12.21 \'2.'6\ 12.41 and 12.51 make A'/'-''(p) a rational function of p such that 
M^*'(0-) = -oo and M'''(1) = (1); so M''^(p) can not be constant on any interval [pi,p2] C (0, 1] and the strict 
behaviour follows. 

Commentary on Lemma 13.21 : 

At this point, it is worth observing that the function a{p) =^ M^^\p) is non-decreasing; strictly positive, 
when p > p^; null, when p = pc and infinetly negative, when p < pc. Again, the reference is Q]. 
Lemma [3.21 above yields 

Corollary 3.3 For each z e N, M^*-* has only one real root in (0,1], denoted by pi^\ in accordace with 
Section \2.4\ 

Lemma 3.4 The sequence of functions | Af*^'^ (p)}^^^, is non increasing, 

ie i < j; i,j eN ^ M'-'\p) > M''^\p), Vp e (0, 1]. In particular a{p) < A//(*)(p);Vf G N, Vp G (0, 1]. 
Idea of Proof: 

The non- increasing behaviour (in i) of |i\f^'''(p)| follows from the First Kind of Coupling ( inequality El 
above) and again from Lemma l3.ll 

□ 

Lemma 13.41 yields 

Corollary 3.5 The numerical sequence \pi^^\ is non- decreasing. Hence limi^oo Pc*"* is well defined 

I J ien 

and p^c^ y i limi^oo Pc"^ ■ Moreover pc pi°"^ > pi^'' , Vi G N 
and 

Corollary 3.6 pi""^ > limi^acPc^- 

Now we turn our attention to the reverse (and more difficult) inequality, viz.pc < limi^oo Pc*"* ■ For that 
consider the following probability spaces: 

• (r2,P) : the abstract probability space, where the iid rvs were defined; 

• (5W,7rW) : the finite probability space of the IMC Y^''\ endowed with its stationary measure tt^*); 

• (f7 X 5(*),P X TT^^') : the product space. 
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In the product space, we define the stochastic processes : 

• I pX^n I , by ;X^^ ((c, a)) {u), V(^.a) G n x S^^) 

. { fx';:^} by -'-'A^f ((c.,a)) 1^'' ^^Xi^^c.), V(a;.a) e il x 5^, 



dcf I ° 

where ryo- G {0, 1}^° is such that rja{2m) = < (t(— 2m) 
starts on state cr e 5'^'^. 



m > 

-1 < m < ; so that the IMC fY^^ 
m < —i 



Now we state 



Lemma 3.7 E 



> M^'\p), yn,ie N, pe (0,1] 



Proof: 

Inequality El implies that V{Lo,a) G x 5<^', pA^i'' ((i^, a)) = > ^■-J('^\uj) = ^^'''^^i'' ((a;, cr)). 

Integrating ( with respect to P x vr*^'-') both sides of the inequality above, yields: 



E 



where the last equality above comes from the fact that the process | 'a'^''} has stationary increments of 
mean 



□ 



Figure 121 below refers to the next lemma 




Figure 3: No reinfection can take place in the region inside the dashed parallelogram. The processes „xi^^ and 
^jj^(oo) indistinguishable there. 



Lemma 3.8 limj^oo E pX^^ = E pX„ for all (fixed) n eN and p e (0, 1]. 
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Proof : 

As the jumps to the right are bounded by +1, it follows that Vcj G Q., pX^'' < n => p^n' — m < n ~ m. 
However, for a site {n,m) to be infected by force, we must have pX'^^ — m > 2i fPef. ITTT b) above). Hence, 
if n — m < 2i, then pX„ — m < 2i,\lui € and site (n,m) will not be infected by force, for any w € . 
Accordingly, it can be proved by induction on n that " n — m 

sites {n,m) such that n~ m <2i have the same infection state regarding the processes pXi°°' and pXl'-*. ( See 
FigureOabove for the case; i = 4.) 
Hence, for any (fixed) time n, 

pXi°°' < pXi'' =^ pXi°°' (m) = 0, Vm > n - 2i (13) 
The reasoning behind ( 1131 is as follows: 



p^n > n — 2i pX^ ^ p^n j — p^n' p^n j ^ "'^ 

pXi'' =^ sup{m G Z„ : pXi'^(m) = 1} = sup{m G Z„ : pXi°°'(m) = 1} pXi°°' =4> 
{ pX« > n - 2i} C { pX« = pXi°°'} ^ { pX« / pXi°°'} C { pX« < n - 2i} C { pXi°°' < n - 2i} 

( Unil follows identifying | pXi'' / pTi°°' | with | pXi'^ > pXi"' | and | pTi°°' < n - 2i| with 

I pXi°°' (m) = 0, Vm > - 2i| . 

Now, observe that, if the infection is not present on the set {(n, m) G Vn : m > n — 2i}, all the paths joining 
its sites to slice Vo must be obstructed somewhere. In particular all the (i+1) straight lines joining site (0, —2j) 
to site {n,n — 2j), < j < i must be interupted at some point (Figure above). As these lines are made of 
different, independent bonds, the probability of this event equals (1 — p^Y'^^ , and we have 

P pXi"' / pXi"') < P pXi°°>(m) = 0, Vm > n - 2i) < (1 - p")'+' 
So that pX^*'' — ^ pXi°°' in probability, and there is a sub-sequence (ifc)fcgN such that pX^''^ — *U pXi°°''a s 



11] Theorem 7.6] . As the whole sequence | pX^*' | is non increasing (in i) , we must have pX^*' — ^ pXi°°' 
well, and we can apply the Monotone Convergence Theorem to conclude that E |^ pX^^'j \' E |^ pX^°°' j 



a s as 



Vn G N. 
□ 

Now we can prove the Convergence Theorem: 
Theorem 3.9 pi*^ Pc 

Proof : Suppose that limi^oo pi*' < Pc, so that we can choose p such that limi^oo pi*' < p < Pc- Then Lemma [3. 2 1 
and Corollary 031 ensure that Vi G N, M'^^[p) > M^'\pi^^) and Lemma ITtI that E [ pXi"'j > nM^'^ (p) > 
0,Vi G N,Vn G N. Thus 

lim E [ pXi''l > 0, Vn G N (14) 

(inequality |S| ensures that this limit is well defined) 
By the other side. 



p < Pc — — > —oo a s Lemma 13. II and commentary on Lemma l3.2l 

3n:E[pX^] <0 



£ I p-^n I SUl^ _QQ Fatou's Lemma 
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lim E [pX^] = E[pX^] < 



Applying Lemma 13.81 for n yields 

lim E I pXi-l = E 

which contradicts inequality 1141 above. Hence, limi^oopi*' > Pc and the theorem follows from Corollary 13.61 

□ 

4 Numerical Calculations 

4.1 Algebraical Determination of the Critical Probabilities 

Theoreni l3.9l ( The Convergence Theorem) is of theoretical interest by itself. However weaker results such 
as Corollaries 13.51 and 13.61 already indicate that each pc^ , i S N, is an improved lower bound (regarding 
its predecessors) for the critical probability of Oriented Percolation. 

Although tacitly present in Sections l2.3.3l and l2.4l above. we present below the algorithm for calculating 
the critical probabilities pc \ i £ N exactly. 

According to Definitions 12.21 12.31 12.41 and 12.61 the critical probabilities p^c \ i G N, may be 
determined in the following steps: 

(i) Determination of the Probabilities P(/^„j ^.^ (Def. 12.3(1 in terms of q^= 1 — p : 

(i) 

elementary combinatorics show that the probabilities PA , , < m < 2' — 1, A; G N, may be 



P^mM) = ■ Polinomial in q (15) 

(ii) Determination of the Transition Matrix (q^'A of the IMC „Yj 



expressed as polinomial functions of q, ie 

{l,m.k) {l,m.k) 

the transition probabilities defined in ( lll|) may be expressed as 

oo 

9^2 = E^^U)('^) (16) 

fc=0 

As the numerical sequence ^P('"/|„ ^ geometric progression (ratio g^) starting from the 

(i + 2)*^* term, the transition probabillities q^^^ may be expressed as rational functions of q. As a 
matter of fact, these probabilities are always polinomials in q: 

(111 = 9iliq) ■■ polinomial in q 

(iii) Determination of the Stationary Measure (vr*) of the IMC pY^^ : 

the combinatorial calculus that lead to equation ( I15|l above show that the transition probabilities 
g[^(q),0 < l,m < 2* — 1, are strictly positive for q e (0, 1), so that the transition matrix ^g^^^ 
is irreducible and aperiodic. Hence tt* is a well defined probability measure on 5* and can be 
algebraically determined from ^<Z;^^j^ ; resulting that 

'nf'{q) : a strictly positive rational function of q, q £ (0, 1), < / < 2* — 1 
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(iv) Determination of the Mean Jump of the SDTCP pX^^ on State ai : 

Definition 12.41 and tlie fact tliat tlic numerical sequence ( -P;ll''t('J') ) is essentially a geometric 
progression yield that 

Af;''*^ (q) : a rational function of q, < / < 2* — 1 

(v) Determination of M^'\ the Mean Jump of the SDTCP pX^^ : 

Definition 12.51 and steps (iii) and (iv) above ensure that 

M^^\q) : a rational function of q 

(vi) Determination of the Critical Probability pi'^ : 

according to CoroUarv 13.31 and Definition 12.61 pc*'' is the only real root of M^^\p) on the interval 
(0, 1); so that it is the only real root of the polinomial in the numerator of M^^^{p). Its calculation 
thus can be done numerically. 

In the sequel, we employ the algorithm described above for calculating the first critical probability in 
algebraical terms: 

4.1.1 The Critical Probability of Zero*'* Order , : 

In this case — 1, so pY^^^ is trivial: 

oo 

M(0)(<z) - ACiq) = ^(1 - 2k)4%,(q) = l.p-l.q{l - q') - 3.qHl - q") - 5.g^^(l - q") - . . . 

k=0 

= (1 - g) - -z(l - q^). [1 + 39^ + 5q^ + Tg*^ + . . .1 = (1 - g) - qil - q^)- 



fe=0 

(l-g2)2 

q + q^ l-2q- q^ 



1 - g2 1 _ g2 

Hence, 



M(°)(g) = 0^1-2g-g2^0^ p^^^ =2-^2^ 0.58579 
4.1.2 The Critical Probability of First Order, p^c^ : 



q-q^ + q^ I - q + q^ - q"^ 
g2 1 — g2 



(1). ^ {q^,l-q + q^~q'^'' 
T^ l W - 



1 - g + g2 + q3 _ 



1-q' ' ' 1-q^ 

1 - 3g + 2g2 - eg-* + + iq^ 



1 - g + 2g3 - 2^4 - g5 . 
Hence, 

M^^\q) = ^ p(^) = 0.604233 . 
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4.1.3 The Critical Probability of Second Order, p. 



(2) 



-q'+q"- 

q* 



+ q 



q-2q^+ q-' + q'' - q" 

q-q' + q'' -q" 

2q^ - 3q^ + + 9^ - 9** 
2 4 

9-9 



1 - 29 + 9^ + 9* - 9° 
9 - 29^ + 29^ - 9'' + 9^^ - 
9 - 29^ + 29^ -9*^+9^-9'^ 
9^-9^ 



q' - 29* + 9' 
1 - 29 + 29^ - 9^ - q-' 
-q-2q^ + 3q^ - q'" - 
1 - 29^ + q" 



'(9) 



'(9) = 



-29* + 29^^ - iq" + 99^ - 149** + 159° - 99^° + 29"- 



(9) 



4'' (9) 



-1 + 39 


-69^ 


+ 893 - 


179^ + 309-'> - 4498 + 469^ 


- 2O9** - 


1798 + 3891" - 


329I1 


+ I39I2 


- 29I3 






-9=^+29 


- 29* + 9° - 49" + 99^ - 


59** - eq" + 12q^° - 8q" 


+ 29!^ 




-1 + 39 


- 692 


+ 893 - 


179^ + 309-'i - 4496 + 469^ 


~ 2O98 - 


179O + 389I0 - 


329I1 


+ I39I2 


- 29I3 








-9^ + 9^ + 29^ - 99*^ + 159^ - 149* 


+ 89^ - 291" 








-1 + 39 


-692 


+ 893 - 


179'' + 309=5 - 4496 + 469'' 


- 20gS - 


179*' + 389I" - 


329I1 


+ I39I2 


- 29I3 


-1 + 39 


-49= 


+ 59^ - 


139'' + 259^ - 279'' + 139^ 


+ 139** - 


349" + 379^° - 


269^^ 


+ 119' = 


-29" 


-1 + 39 


- 692 


+ 893 - 


179^' + 3O95 - 44g<5 + 469^ 


- 2O9S - 


1799 + 389I" - 


329II 


+ I39I2 


- 29I3 



M<^'(9) ■ 



Hence, 



My 



1 - 29 - 59 + 49 , 1 - 29 - 39 + 29 , 1 - 29 - 9 - 29 + 29 



-1 + 59 - 119^ + 169^ - 259* + 529^ - 759*^ + 969^ - 589** - 699^ + 1529^" - 1119^ 



1 - 29 - 9-j 
- 59^^ + 749I 



■ 499" + 109^'' 



(1 - 92) (-1 + 39 - 692 + 893 - 179* + 3O95 - 4496 + 469^ - 2898 - 1799 + 389IO - 329I1 + 1312 - 29I3 



M^^\q) = : 



0.614187 



4.2 Numerical Determination of the Critical Probabilities: 

The combinatorial calculations leading to the transition matrix and to the mean-drift vector 

M(z); described in Section 14. II above, yet fastidious for humans, are taylor-made for computers: while it 
took us a whole afternoom for calculating algebraically (16 polinomial entries), a FORTRAN 77 

program, running at 600 MHz, calculated ('Z*-^'');^ (4^ polinomial entries) in less than one minute time! 

The greatest problem in writing a computer program version for the algorithm described in Section 
I4.1l arises precisely in step(iii), viz. the determination of the stationary measure tt^'^ (in algebraical terms), 
starting from the transition matrix (9*-*'');^ (l)- True, it can be done straightforwardly, acting on the 
polinomials as if they were numbers. However, if we don't take into account the eventual and probable 
simplifications that may take place, the degree of the outcoming polinomials soon get unmanageble. And 
finding out these simplifications is far from obvious. 

Therefore we have adopted a different approach: 

(i) After having algebraically calculated the transition matrix (<?'•*•') (g) and the mean-drift vector 

Ml'\q), the program calculated a numerical transition matrix and a numerical mean-drift vector 
for a decreasing sequence of numerical values of q. 

(ii) From each numerical transition matrix, a numerical stationary measure was obtained by solving a 

system of linear equations. It is worth-while mentioning that partial pivoting was of fundamental 
importance in keeping numerical errors under controU. 



(iii) The numerical 



-drift was then obtained, performing the inner-product of Definition 12.51 



(iv) According to Definition 12.61 ( 1 — Pc ) lies between the last value of q for which M^^\q) is 
negative and the first value of q for which Af'^*^(g) is positive. 
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This approach has the onus of introducing numerical erros (basically in step ii above), from which we 
hitherto have got rid of. So as to produce reliable numerical data for the critical probabilities, rigorous 
upper bounds for the numerical errors must be provided. In doing so, we followed the set-up of foward 
analysis described in which consists fundamentally in providing partial upper bounds after each 
arithmetic operation performed in the flow of numerical calculations leading to the final numerical result. 
This rather crude approach has the advantage of providing a secure upper bound for numerical errors, 
regardless of particular features the linear system may present, as it is often the case in the so called 
backward analysis ( 

In the sequel, we present the numerical data produced following steps i-iv above^ : 



i=l 



p 




AM(i)(p) 


0.604231 
0.604232 
0.604233 
0.604234 


-9.15E-6 
-5.27E-6 
-1.40E-6 
+2.47E-6 


6.2E-16 
6.2E-16 
6.2E-16 
6.2E-16 



1=2 



1=3 



p 




AM(2) (p) 




P 






0.614185 


-l.OlE-5 


6.1E-15 




0.620203 


-1.36E-5 


2.5E-14 


0.614186 


-5.85E-6 


6.2E-15 




0.620204 


-8.96E-6 


2.5E-14 


0.614187 


-1.57E-6 


6.1E-15 




0.620205 


-4.32E-6 


2.6E-14 


0.614188 


+2.71E-6 


6.2E-15 




0.620206 


+3.15E-7 


2.6E-14 



1=4 



p 




AM(4)(p) 


0.624210 
0.624211 
0.624212 


-9.52E-6 
-4.57E-6 
+3.71E-7 


8.0E-14 
8.0E-14 
8.0E-14 



i=5 



i=6 



p 




AM(5)(p) 




P 




AM(«)(p) 


0.627064 


-1.31E-5 


1.9E-13 




0.629201 


-1.55E-5 


3.8E-13 


0.627065 


-7.85E-6 


1.9E-13 




0.629202 


-9.98E-6 


3.8E-13 


0.627066 


-2.63E-6 


1.9E-13 




0.629203 


-4.51E-6 


3.9E-13 


0.627067 


+2.60E-6 


1.9E-13 




0.629204 


+9.66E-7 


3.8E-13 



1=7 



i=9 



p 




am(^) (p) 


0.630863 
0.630864 
0.630865 


-l.OlE-5 
-4.41E-6 
+1.30E-6 


7.2E-13 
7.2E-13 
7.2E-13 



P 




AM(**)(p) 


0.632192 
0.632193 
0.632194 


-8.53E-6 
-2.60E-6 
+3.32E-6 


1.3E-12 
1.2E-12 
1.3E-12 



P 




AM(^>(p) 


0.63326 
0.63327 
0.63328 
0.63329 


-1.29E-4 
-6.79E-5 
-6.68E-6 
+5.45E-5 


2.0E-12 
2.0E-12 
2.0E-12 
2.0E-12 



Based upon the data displayed above, the first ten critical probabilities may be determined: 



i 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Pc' 


0.585787 


0.604233 


0.614187 


0.620205 


0.624211 


0.627066 


0.629203 


0.630864 


0.632193 


0.63328 



So that we can establish the inequality : 



Pc > 0.63328 



"the calculations were carried out on double precision. 



'^|M^*^(p) — -^^^Tn(?')l — ^M^^Hp) '■) where -^^^^(p) denotes the numerically calculated mean drift of i^^ order. 
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5 Simulations 



The boolean and inductive features of definition I2.1l make it extremely suitable for computer simulations: 
at each step a FORTRAN77- written computer program determmines the value of pX^^ /n, the right edge 
mean speed at moment n. In accordance with lemma 13. II this sequence converges (a.s. in n ) to MW(p), 
the right edge mean speed. Therefore running the program for different values of p, one can estimate the 
value of pc^ observing the height of the plateau established with the increase of n, according with the 
rule: 

(a) height of plateau < ^ p < pi^^ , 

(b) height of plateau > =^ p > p^f' . 

The two figures displayed below ilustrate the use of this technique when i = 1000: 



Xn/n 

i = 1 ooo 




log n 



Figure 4: Independent (not coupled) trajectories of the processes .64450^n ^ .64452^n /n. 



Xn/n 

i=i ooo 




log n 



Figure 5: Enlarged detail of the figure above . 
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Table ^ below sumarizes the numerical results obtained by means of the technique for different values 
of i : 



i 


Pc 


5 


0,627 


6 


0,629 


9 


0,6332 


20 


0,638 


40 


0,641 


100 


0,643 


200 


0,6438 


1000 


0,64451 



Table 1: Critical probabilities obtained by means of simulation. 

It is interesting to observe that the plateau pattern exhibited in figure ^is present even in the vicinity 
of criticality, ie p ~ pi*' and i oo. Therefore the critical probabilities pc*"* (even for large values of i 
and so very close to criticality) can be estimated within the desired precision simply by increasing the 
value of n, ie running the program until the plateau pattern becomes clear. We belive that this feature 
distinguishes our technique from the usual Monte Carlo simulation methods, wherein instability appears 
near criticality and the use of scaling techniques is called for. 
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